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We present a new method for analysing directed flow, based on a three-particle az- 
imuthal correlation. It is less biased by nonflow correlations than two-particle methods, 
and requires less statistics than four-particle methods. It is illustrated on NA49 data. 

O ■ 1. INTRODUCTION 
00 

The measurement of the successive harmonics of azimuthal correlations in a heavy 
ion collision |T[ , v n = (e m ^~^ R ^ ^> , where $ r is the impact parameter direction, is of 
paramount importance, since it yields information on the medium created in the collision. 
In particular, a novel behaviour has been predicted for the first harmonic the so-called 
directed flow, at ultrarelativistic energies [[| . 

However, the analysis of directed flow at these energies is highly nontrivial, because v\ is 
very small. Thus, methods based on an analysis of two-particle azimuthal correlations 
are likely to be biased by "nonflow" two-particle correlations: quantum (HBT) correlations 
between identical particles and correlations due to global momentum conservation have 
been shown to be important at SPS Q, while correlations due to minijets may be large at 
RHIC U . On the other hand, methods relying on four-particle correlations [H] , which are 
free from this bias, are plagued by a lack of statistics due to the smallness of V\, although 
they give good results for the analysis of elliptic flow v 2 (see Sec. Q). To remedy these 
shortcomings, we proposed a new method of v\ analysis @, based on the measurement of 
a mixed three-particle correlation, which involves both v\ and v 2 : 

,i(<t>a + <t>b-2<t>c)\ ^ \2 



e H0 o +0 b -20c^ ^ {vxfVi, (1) 

where (f> a , fa, and C denote the azimuths of three particles belonging to the same event, 
and the average runs over triplets of particles emitted in the collision, and over events. 
Once v 2 has been obtained from a separate analysis, this equation yields (t>i) 2 i>2, thus V\. 

Here, we apply this method to NA49 data on Pb-Pb collisions at 158 AGeV. Results 
obtained using the "standard" flow analysis are given in Ref. ||. In our method, the first 
step in the analysis is the measurement of a reference v 2 - Then, an equation analogous 
to Eq. [I] yields an estimate of the integrated Vi, i.e., its average value over some phase 
space region. Finally, restricting <pi in Eq. [I] to a small (p?, y) bin allows one to obtain 
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more detailed, differential measurements of v\ as a function of transverse momentum or 
rapidity. 

2. ELLIPTIC FLOW FROM 2, 4, 6, 8-PARTICLE CORRELATIONS 

As stated in the introduction, our method of analysis of directed flow V\ requires the 
preliminary knowledge of an estimate of the elliptic flow v 2 , integrated over some phase 
space region. Of course, this estimate must be obtained by analysing the same sample of 
events from which one wants to extract V\. 

In practice, the average over phase space is a weighted average: 



where w 2 is the chosen weight. In order to reduce statistical fluctuations, w 2 must be 
larger for particles with stronger elliptic flow. 

The value of (w 2 v 2 ) is obtained using the cumulant method described in Ref. ||: one can 
extract estimates of (w 2 v 2 ) from cumulants of multi- (2-, 4-, 6-. . . ) particle correlations. 
While two-particle methods are equivalent to the standard flow analysis, higher orders 
are essentially free from nonflow effects. They were first used in analysing data obtained 
by the STAR Collaboration at RHIC 0. 

In Fig. |I|, we present application of the method to NA49 data, and show the dimen- 
sionless quantity 



as a function of centrality for charged particles, where we have used w 2 = pr- We display 
estimates using cumulants of two-, four-, six-, and even eight-particle correlations [|], 
as well as the corresponding quantity for charged pions obtained from the "standard" 
subevent method 0. 

It is quite remarkable that the four-, six-, and eight-particle estimates all agree: this 
supports the idea that they are indeed free from nonflow effects, and correspond to a 
genuine collective motion in the direction of the impact parameter. Moreover, these 
multiparticle estimates show a slight discrepancy with the two-particle values, as expected 
if nonflow correlations are sizable ||. Please note that the statistical uncertainties on 
high order cumulants remain reasonably small, especially for midcentral collisions. In 
the following, our reference elliptic flow value will preferably be the estimate from the 
four-particle cumulant, which is a priori the most reliable since it is free from nonflow 
correlations and has a smaller statistical error than higher order estimates. 



3. INTEGRATED DIRECTED FLOW FROM 2, 3, 4-PARTICLE CORRE- 
LATIONS 



The next step in the analysis is to determine the average value of directed flow over 
some phase space region. As in the case of elliptic flow, we perform a weighted average 




(2) 



(w 2 v 2 ) 



(3) 





(4) 
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where stronger weight is given to particles with stronger directed flow. In the NA49 
analysis, we used a rapidity dependent weight w\ = y — ycM, where ycM is the centre-of- 
mass rapidity. 

This weighted average is obtained from the following three-particle correlation: 

(w 1 (a)w 1 (b)w 2 (c)e i{ ^ b - 2 ^ ) ) = (w 1 v 1 ) 2 (w 2 v 2 ) , (5) 

where w 2 is the same as in Sec. |2|. Using the value of (w 2 v 2 ) obtained in Sec. 0, we thus 
derive (wiVi), up to a global sign. 

In practice, the left-hand side of Eq. || is constructed using a generating function for- 
malism detailed in Ref. [0]. This procedure is a very efficient way to sum over all pos- 
sible triplets of particles, and also to remove automatically the effects of slight detector 
anisotropics. 

The method was applied to NA49 data. In Fig. |2|, we present as a function of centrality 
the dimensionless quantity v 1 = (wiVi) / J ((u>i) 2 ), analogous to Eq. |, together with two- 
and four-particle estimates obtained with the method of Ref. [|]. It is worth noting that 
the statistical uncertainty on the three-particle estimate is barely larger than that on the 
two-particle value, while the systematic error due to nonflow correlations (not included 
in the plot) is a priori much smaller. 
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Figure 1. Weighted elliptic flow v 2 of 
charged particles function of central- 
ity (central left, peripheral right) for Pb-Pb 
collisions at 158 AGeV. 



Figure 2. Weighted directed flow v\ of 
charged particles function of centrality 
(central left, peripheral right). 



4. DIFFERENTIAL DIRECTED FLOW 

Finally, we want to obtain detailed measurements of vi, as a function of px or y and 
for each particle type. To derive vi(pT,y), we use a three-particle correlation analogous 
to Eq. H where the particle labeled a is restricted to a narrow (pr,y) bin: 

{ Wl {b)w 2 {c)e^ + ^- 2 ^) = (w lVl ) (w 2 v 2 ) Vl {p Ta ,y a ). (6) 

With the previously derived (w 2 v 2 ) and (w\Vi), this equation yields t>i(pr, y)- We illustrate 
this method on the differential flow of charged pions in midcentral collisions, shown in 
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Figs. |3] and f| together with two two-particle estimates, either uncorrected or corrected 
for the effect of momentum conservation [|1(J . 



At high px, V\ from three-particle correlations is consistent with the two-particle value 
corrected for momentum conservation, but significantly lower than the uncorrected one. 
This shows that correlations from momentum conservation, which are large, are automati- 
cally removed in our method. This is also reflected by the behaviour at midrapidity, where 
V\ vanishes, as it should by symmetry, while the uncorrected two-particle estimate does 
not. The three-particle estimate also vanishes more smoothly at pt = than two-particle 
estimates which may be biased by HBT correlations [1. 
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Figure 3. Directed flow v\ of charged pions 
as a function of transverse momentum for 
midcentral Pb-Pb collisions at 158 AGeV. 



Figure 4. Directed flow v\ of charged pi- 
ons as a function of y — i/cm for midcentral 
collisions. Open points are reflected with 
respect to midrapidity. 
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